Measures of entanglement in multipartite bound entangled states 
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Bound entangled states are states that are entangled but from which no entanglement can be 
distilled if all parties are allowed only local operations and classical communication. However, in 
creating these states one needs nonzero entanglement resources to start with. Here, the entanglement 
of two distinct multipartite bound entangled states is determined analytically in terms of a geometric 
measure of entanglement and a related quantity. The results are compared with those for the 
negativity and the relative entropy of entanglement. 
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Introduction — We are motivated to study the quantifica- 
tion of entanglement for the basic reason that entangle- 
ment has been identified as a resource central to much of 
quantum information processing (see, e.g., Ref. 0]). To 
date, progress in the quantification of entanglement for 
mixed states has resided primarily in the domain of bi- 
partite systems ^ . For multipartite systems in pure and 
mixed states, the characterization and quantification of 
entanglement present even greater challenges. 

Among all entangled states, there is a peculiar class 
of states, bound entangled states, originally discovered 
in the bipartite setting, where some connection to zero 
negativity was identified Q. These are states that are 
entangled, but from which no pure entangled state can 
be distilled if all parties are allowed only local operations 
and classical communication (LOCC) . The distillable en- 
tanglement (-Ed) is thus zero. Bound entangled states 
can be either bipartite or multipartite, the latter possi- 
bly exhibiting more structure than the former. However, 
it does take nonzero entanglement to create bound en- 
tangled states under LOCC. It is thus desirable to see 
how much entanglement is needed. 

However, the two most important measures — 
entanglement of distillation (-Ed) and of formation 
{E-p) — have so far been limited to bipartite settings, as 
there are ambiguities in generalizing them to multipar- 
tite systems 0] ■ In order to explore multipartite settings, 
it is thus, on the one hand, necessary to lay down bounds 
on the entanglement content for distillation and forma- 
tion. On the other hand, applying other measures, such 
as the relative entropy of entanglement (Er), to mul- 
tipartite states may prove helpful in quantifying entan- 
glement. Recently, a multipartite entanglement measure 
based on the geometry of Hilbert space has been pro- 
posed 0, IE > and has been applied to several bipartite 
and multipartite cases. The merit of this measure is that 
it is suitable for any-partite systems with any dimension, 
although determining it analytically for generic states re- 
mains a challenge. 

In the present paper, we study the entanglement con- 



tent of two distinct bound entangled states: Smolin's 
four-party unlockable bound entangled state 0, and 
Diir's iV-partv Bell-inequality-violating bound entangled 
states [ill. For each, we determine analytically their ge- 
ometric measure of entanglement and a related quantity. 
Under certain circumstances, these give lower bounds on 
their multipartite E-p- In addition, we make conjectures 
concerning their relative entropies of entanglement. Al- 
though quantities such as the geometric measure or the 
relative entropy of entanglement may not be able to re- 
veal the exact nature of bound entanglement, they nev- 
ertheless quantify for these bound entangled states the 
content of entanglement that is unextractable. 

Geometric measure of entanglement — We begin by briefiy 
reviewing its formulation. Consider a general n-partite 
pure state (expanded in the local bases {cp}})'- 
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As shown in Ref. f^, its closest separable (i.e., product) 
pure state (with i being the party index) 

1^)^ ® ® (X c«|e«)), (2) 

1=1 1=1 V-"^ — 'Pi J 

satisfies the condition (and its complex conjugate) 



p\---Pi---P' 



^P\P^---PnP\ ^Pi ^Pn ~ ^Pi ' 



where the eigenvalue A e [—1,1] is associated with 
the Lagrange multiplier enforcing (010) = 1, and the 
caret denotes exclusion. Moreover, A is the cosine of 
the angle between \ip) and |0); the largest one Amax, 
which we call the entanglement eigenvalue, corresponds 
to the closest separable state, and is the maximal over- 
lap: Amax(IV')) = niax0 K^IV') I, where \(f>) is separable 
but otherwise arbitrary. E^in^ = 1 — A^-.^xd'/')) "^as de- 
fined to be the geometric measure of entanglement 
for state {ip), and it measures the degree of inseparability 
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via the squared sine of the angle away from the closest 
separable pure state. 

The extension to mixed states can be built upon the 
pure-state theory and is made via the use of the convex 
hull construction, as was done for Ep Q . The essence is a 
minimization over all decompositions p = J2iPi |V'i)('0i| 
into pure states: 



>.{p) = min } 



(4) 



This convex hull construction ensures that the measure 
gives zero for unentangled states; however, it also compli- 
cates the task of determining mixed-state entanglement. 
Egi^2 was shown to be an entanglement monotone 0] 
(i.e., the measure does not increase under LOCC), hence 
is a good measure of entanglement. As there is no explicit 
generalization of Ep to multipartite states, we shall cal- 
culate i?sin2 analytically for two bound entangled states, 
Smolin's and Diir's. Because E-p{p) is the minimum av- 
erage ehit to create a single copy of p, we can regard 
_Esijj2 (p) as the minimum average degree of pure insepa- 
rability needed to realize the state p. 

In bipartite settings, it is known (llj that Eyi{p) < 
Ef{p), and that for pure states Itp), -Er(V') = Ep{ip). 
It is also known |l2l| that for any bi- and multipartite 
pure state Eio^^i.'^) = -2 loga Aniax(V') < -Br('0). 
Together with the inequality (1 — x^) logj e < — 2 log2 x 
(for < X < 1), one has 



Vi E^in^- (-04 ) l0g2 e<^p^ flog, (ll^i) <^,Pt Ef{jP, 



and thus (logg e)£'sin2 (p) < E\og^{p) < Ep{p) for any bi- 
partite state p. If the generalization of Ep to multipartite 
states maintains the property that Ep{il}) > E-[i{ip) then 
the inequality (log2 e)-E'si„2 < Eiog^ < Ep will continue to 
hold for multipartite mixed states. We remark that 



- 2l0g2 Amax(V'j)) (5) 

is not an entanglement monotone |l2l |. However, we see 
that both (log2 e)i?si,j2 and Sio^^ could serve as lower 
bounds on multipartite entanglement of formation. 

We now turn to the calculations of entanglement for 
the two bound entangled states: Smolin's and Diir's. 

Smolin's four-party unlockable bound entangled state — 
This is a four-qubit mixed state 



ABCD 



iE(l*^)(*»l)AB«(l*»)(*»l)cD' (6) 

1=0 



where the |^')'s are the four Bell states (|00) ± |ll))/\/2 
and (|01) ± |10))/y2. Now, the state p^^CD 
conveniently rewritten as 



ABCD 



1 ^ 

-T 

j=0 



(7) 



where the |X)'s are the four orthogonal Greenberger- 
Horne-Zeilinger-(GHZ) like states: 

|Xo)^-i=(|0000)-f 11111)), |Xi)^i=(|0011) + |1100)), 

|X2)^i=(|oioi)-Mioio)), 1x3) ^-^(10110) + 11001)). 

The most general decomposition of a mixed state p 
into pure states can be expressed as 

M n 

P = '^\V'k){<fk\ with \ifk) = y^^fciVX l^z), (8) 



k=l 



1=1 



where M is an integer not smaller than n, the number of 
orthonormal eigenvectors {|^i)} (with nonzero eigenval- 
ues {Xi}) of p, the |(^)'s are unnormalized , anadU satisfies 



Thus, the most general pure state 



that appears in the decomposition of Smolin's state is 



\Vk) 
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E 9 ^fc^l^i)- 



(9) 



i=0 



Our goal is to minimize X^fcP*: ^puio(l'/'/c)) over all pos- 
sible Zi's, where iSpuro is some pure-state entanglement 
(i?si„2 or fiog2 in our considerations), pu = (<^fe|<^fc), and 
|(/9fe) is the normalized state \Lpk) = \(pk)/\/Pk- Making 
a general minimization for an arbitrary mixed state is 
extremely difficult. However, for the mixed state p^^'^^ 
we shall show that the decomposition in Eq. iQ does in- 
deed minimize the average entanglement over pure-state 
decompositions. As in Eq. ©, It^s) can be explicitly writ- 
ten as |</?) — X]i=o l'''^*)' "^J^cre the g's are non- 
negative, satisfying ^^qi = 1, and the (/)'s are phases. 
For fixed q's, the state has a maximal entanglement eigen- 
value when all phases are zero. We shall show shortly 
that its maximal entanglement eigenvalue is l/\/2, which 
is achieved by the |A")'s. 

The entanglement eigenvalue of the state = 
Y^i=o V'i'ilXi) is the maximal overlap with the separa- 
ble state 1$) = (^j=i{ci\0) -hsi\l)), where a = cos 9i and 
Si = sin 61^ with < Oi < n/2. Thus 



{^\ip) = VW2 (C1C2C3C4 -f S1S2S3S4) -I- y/qi/2 (C1C2S3S4 
+S1S2C3C4) + V 172/2 (C1S2C3S4 -I- S1C2S3C4) 
+ \Jqzl2 (C1S2S3C4 + S1C2C3S4), 

which has maximum I/V2. To see this, use 
the Cauchy-Schwarz inequality, treating as one vec- 
tor { v^go/S, \Jq\l2, ^Jqil2, \Jqzli) (whose modulus is 
1/a/2 ) , and the corresponding coefhcients as another vec- 
tor (whose modulus can be shown to be no greater than 1; 
see App. A). The states \Xij clearly saturate this bound; 
hence: 



(10) 
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This suggests that although bound entangled, Smolin's 
state has a very high degree of entanglement, the same as 
that of a 4-partite GHZ state. This high degree of entan- 
glement seems to manifest in some bipartite partitioning, 
e.g., {A:BCD} (as we discuss below). 

We conjecture (and later prove) that Eji = 1 for this 
state and one of its closest separable mixed states is 

^(|0000)(0000| + |1111)(1111| + |0011)(0011| + |1100)(1100| 
8 

+ |0101)(0101| + |1010)(1010| + |0110)(0110| + |1001)(1001|). 

in 

We remark that the negativity M 2] (a value used to 
quantify the degree of bipartite inseparability of states 
and defined as twice the absolute sum of negative eigen- 
values of the partial transpose (PT) of the density matrix 
with respect to some bipartite partitioning) is zero for 
any 2/2 partitioning, e.g., {AB : CD}, but nonzero for 
1/3 partitioning, e.g.,{A:BCD}. Specifically, A/aiBCD — 
1 but 7Vab:CD = 0. 

Let us now turn to Diir's bound entangled states. 
Diir's N -party bound entangled states — Diir ^3 found 
that for > 4 the following state is bound entangled: 

PN ^ (^i*g>(*g| + ^ E (^fc + h)^ , (11) 

where |*g) = (|0®^) + e'°"^\l'^^)) /V2 is a Af-partite 
GHZ state; Pk = \uk){uk\ is a projector onto the 
state \uk) = |0)i|0)2 . . . |l)fc . . . |0)iv; and Pfc = \vk){vk\ 
projects on to \vk) = |l)i|l)2 ■ • • |0)fc . . . |l)Ar. For A^ > 
8 this state violates the Mermin-Klyshko-Bell inequal- 
ity violation was pushed down to A^ > 7 by Kasz- 
likowski et al. ^3 ^ three-setting Bell inequality; it 
was pushed further down to A^ > 6 by Sen et al. for a 
functional Bell inequality. The phase ajv in IvI/g) can be 
eliminated by local unitary transformations, and hence 
we shall take a^v = in the following discussion. 
In fact, if we consider the family of states 

1 - ^ 

Pa,(x)^x|vI'g)(*g| + ^E(^'= + ^0' (12) 

fc=i 

we find that for A > 4 the state is bound entangled if 
< a; < 1/(A^ + 1) and is still entangled but not bound 
entangled if x > 1/(A^ +1)- This can be seen from the 
fact that the negativities of pn{x) with respect to the 
two different partitions (1 : 2 • • • A^) and (12 : 3 • • • A^) are 

AAi,2...Jv(pjv(x)) = max{0, [{N + 1) x-l]/N}, (13a) 
J^12:3-n{pn{x)) ^ X. (13b) 

By applying arguments similar to those used to calcu- 
late entanglement for Smolin's state, we have that the 



general pure state in the decomposition of pn{x) is 

N 
k=l 

where the g's and r's are non-negative and satisfy 
Sfc(9fc + ^fe) = 1- In this family, the state with the least 
entanglement (or maximum Amax) for fixed {y, qk,rk} is 
the one with all phase factors zero: 

N 

|*(y,{(7,r})) = \/y|*G) + v/l^E(\/^l"»)+V^I"*))- 

k=l 

Next, we ask: For fixed y, what is the least entanglement 
that the above state can have? Take a separable state 
of the form |$) = (8i^i(ci|0) -I- s^jl)); its overlap with 
I *(y is then 

(*|$) = ,/jJj2{ci---CN + Si---SN) 

N 

+ \/l-yy^^{y/qkCi • • • Sfc • • -CN + y/r^Si ■■ - Ck-- - Sat). 

fe=l 

This can be shown to be no greater than ■\/(2 — y)/2, 
again by a Cauchy-Schwarz inequality, taking 

[^,y{l-y)qk},W{i-y)rk}} 

as the first (2A^ + l)-component vector (with modulus 
y/ [2 — y)/2) and the corresponding coefficients as the sec- 
ond one (whose modulus can be shown to be no greater 
than 1 for A^ > 4; see App. A). The bound can be satu- 
rated, e.g., by 

\i'±,uAv)) = V^I*g) ± VT^\uk), (14a) 

\^j±.vAy)) = V^I^g) ± V^\vk), (14b) 

for which A,„ax(2/) = V(2-y)/2 [ij. As 1 - Al,^^{y) is 
linear in y and —2 log2 Aniax(y) is convex in y, one gets 

X 2 

Esin^{pN{x)) = -, £iog,(pAr(x)) = log^ ^— ^ , (15) 

and one of the optimal decompositions is 
1 ^ 

^^(^)"4/V^^ ^ IV'a,/3,fe(a;))(V'a,/3,fc(a;)|. (16) 

k—l a— ± P—u,v 

The above calculations show that for pn(x), the entan- 
glement depends on the portion x of the GHZ state in 
states \'ipa,t3,kix)){ipa,f3,kix)\ and it never becomes zero 
unless there is no GHZ mixture. 

We conjecture that, for A^ > 4, pn{x) has Eji{x) — x, 
with one closest separable mixed state being 

, _ w 

|(|0..0)(0..0| + |1..1)(1..1|) + ^ E (^fe + ^fc)' 

k=l 
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which seems plausible as (|0..0)(0..0| + is a 

closest separable mixed state to |^g)- 

Concluding remarks — We have presented analytical re- 
sults on how much entanglement is bound in two distinct 
multipartite bound entangled states. The measure we 
have used to quantify their entanglement is the geomet- 
ric measure of entanglement (GME), whose construction, 
similiar to the entanglement of formation (Ep), is via the 
convex hull. In contrast to GME, Ep has not been explic- 
itly generalized to multipartite states, and hence is still 
unavailable for these bound entangled states. However, 
under the circumstances discussed previously, the results 
for Esin2 as well as a related quantity, Siog^ ' might provide 
lower bounds on E-p. For the Smolin state, its bound en- 
tanglement is as large as that of a four-partite GHZ state, 
whereas that for Diir states is related to the portion of 
the A^-partite GHZ state. For each case, an optimal de- 
composition is given. Furthermore, we have conjectured 
that the relative entropy of entanglement (-Er) for the 
Smolin state is unity (proved below), whereas we conjec- 
ture that E^ for Diir's state is equal to the portion that 
is iV-partite GHZ. 

For Smolin's state we can establish its Ep, E^, E^ , and 
i?si„2 for certain bipartite partitionings. For example, if 
we group the four parties ABCD in two, A:BCD, we can 
write the state as 

with the 3-qubit states of BCD mapped on to the 8-level 
system (000 0,001 1, ...,111 ^ 7), involving the 
locally orthogonal and convertible states (by BCD) 

l^o) = (|00)-I-|17))/V2, = (|03) + |14))/V2, 

1^2) = (|05) + |12))/\/2, 1X3) = (|06) + |ll))/%/2. 

In order to find the entanglement of this bipartite state 
(in ® C*), we need to consider the entanglement of 
the general (properly normalized) pure state 

^ — 

that appears in the pure-state decompositions. In fact, 
regardless of the values of the x^'s, this pure state has a 
reduced density matrix (tracing over BCD) of the form 
(|0)(0| + |1)(1|) /2. This shows that p^-BCD j^^^g ^ 
i?gi„2 = 1/2, and Siog^ = 1- ^^.ct, there is a general 
result due to Horodecki et al. that Ed — Ep for 
mixtures of locally orthogonal bipartite states, e.g., (8) 
C^™ states that are derived from mixing Bell-like states 

|*±) = (|0,fc)±|l,2m^_fc^))/\/2, (18) 

having distinct fc's, where k — 0, l,...,m — 1. As 
< < £^F, we have that Ep^{p^-^^^) = 1 as well. 



What about the original four-partite state pABCDy 
i?R(p^^™) > i?R(p^^^^°), we have Ep^{p^^^^) > 1. 
But we also have that £'r(p^^^'^) < 1, as our previ- 
ous conjecture gives at least an upper bound; we thus 
have that Eb,{p^^™) — 1 and the conjecture is proved. 
Naively, we expect that any arbitrary pABCD j^^g greater 
entanglement than pAiBCD^ However, for the Smolin 
state, they have the same entanglement as quantified by 
both GME and the relative entropy of entanglement. 

Although Diir's bound entangled state violates a Bell 
inequality, it has nonzero negativity under certain parti- 
tionings. One may raise the question: Does there exist a 
bound entangled state that has positive PT (PPT) under 
all partitionings but that still violates a Bell's inequality? 
For example, does an unextendible-product-basis (UPB) 
bound entangled state |0| violate a Bell inequality? We 
shall see shortly that the answer is "no" , at least for the 
three different Bell inequalities 0, 0, Q| mentioned 
earlier. Acfn has shown 18] that if an A^-qubit state vi- 
olates a two-setting Bell inequality then it is distillable 
under certain bipartite partitioning. Using the results of 
Refs. 0, regarding distillability, we can repeat the 
same analysis for the other two inequalities [isl and 
indeed obtain the same conclusion; see App. B. This bi- 
partite distillability then implies a negative PT (NPT) 
under that bipartite partitioning according to Horodecki 
et al. p] Hence, violating these Bell inequalities implies 
NPT under certain bipartite partitioning. Said equiva- 
lently, if an iV-qubit state has PPT under all bipartite 
partitionings then the state never violates these Bell in- 
equalities. This seems to suggest that PPT bound en- 
tangled states are truly bound in nature that cannot give 
deviation from local theories. 
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S. Mukhopadhyay, O. Rudolph and especially W. Diir for 
discussions. This work was supported by NSF Grant No. 
EIAOl-21568 and DOE Grant No. DEFG02-91ER45439. 

Appendix A: In this appendix we sketch proofs of two 
useful inequalities and describe the deriviation of the en- 
tanglement eigenvalue for the states in Eqs. (fTIjl . We 
start with the first sought inequality: 

(C1C2C3C4 + SiS2S3S4)^ + (C1C2S3S4 -I- SlS2C3C4)^ 

+ (ciS2C3S4 -I- SiC2S3S4)^ -|- (C1S2S3C4 + SlC2C3S4)^ < 1. 

We have simplified the notation by using c; = cos 6i and 
Si = sin 9i. By subtracting the left-hand side from 1 
and doing some algebraic manipulation, we arrive at the 
non- negative expression (hence the sought result): 

(C1C2C3S4 - SiS2S3C4)^ -I- (C1C2S3S4 - SiS2C3S4)^ -|- 
(C1S2C3C4 - SlC2S3S4)^ -t- (S1C2C3C4 - Ci 525354)^. 
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The next sought inequaHty is (for TV > 4): 

/at = (ci • • • Cat + Si • • • Sat)^ + 
N 

^ {(Cl • • • Sfc • • -CTv)^ + {Sl ■ ■ ■ Ck ■ ■ ■ sn)"^} < 1. 
k=l 

First, making similar arguments, one can show that 
/4 < 1- One can also show that Jn+i < In- Thus, 
by induction, we have the sought result. 

We now discuss why y/yl'^c) ± -y/l ~ vWk) and 
y/y\^!G)±y^l — y\vk) have as their maximal entanglement 
eigenvalue Amax(2/) = \/(2 — y)/2. As one can make local 
relative phase shifts to transform y^l'i'c) + %/! — y\uk) 
to y/yl^'c) — ^1 — y\uk), they have the same entan- 
glement. The change from ^\'^g) ± -^/l — y\uk) to 
Y^I^'g) ± y/l — is simply a flipping of to 1, and 
vice versa. The mapping from k to k' is just a relabeling 
of parties. Thus, we need only consider the state 

\/^(|oo---o> + |ii---i» + v/T^|io---o). 

As this state is invariant under permutation of all par- 
ties except the first one, and as the coefficients are non- 
negative, in order to find the maximal overlap we can 
make the hypothesis that the closest separable state is of 
the form 

(vp|o> + v^|i>) ® (V^|o> + yT^|i))«^-^ 

We further see that in order for the overlap to be max- 
imal, q must be either 1 or 0. For the former case, we 
can further maximize the overlap to get •y/(2 — y)/2. For 
the latter case, the maximum overlap is ■y/y/2, which is 
less than \/{2 — y)/2 (as < y < 1). Hence, the state 
Y^I^'g) ± VI ~ yWk) has the entanglement eigenvalue 

\/(2-2;)/2. 

Appendix B: In this appendix we analyze the connection 
between violation of three Bell inequalities and bipartite 
distillability as was done in Ref. jl^] for the two-setting 
inequality. It was shown by Diir and Cirac |20l | that an 
arbitrary A^-qubit state p can be locally depolarized into 
the form 

= A+|*+)(vI'+| + Ao-|*o)(*o| 

+ E ^.(i*;>(*;i+i*7)(*7i)' 

while preserving Ag = and A-, = + 

where = (|0»^) ± |l®^))/v4, and the 
|'I'^)'s are GHZ-like states, i.e., the states in Eq. (^HJ, 
unfolded into qubit notation. Normalization gives the 
condition 

Aj + Ao-+2^Aj = l. 
j 



Now define A = A,| — A^, which we assume to be non- 
negative without loss of generality. The condition that 
there is no bipartite distillability for some bipartite par- 
titioning Pj is ^3 

2\j > A. 

Assuming nondistillability for all bipartite splittings, we 
have 

2^A,=l-(A++Ao)>(2^-i-l)A. 
i 

As A(j" + Aq" > a, we have further that 

1 - A > (2^-1 - 1)A. (19) 

For the two-setting Bell inequality considered by 
Aci'n [lil, violation implies A > l/2(^-i)/2. For the 
three-setting Bell inequality considered in violation 
implies A > \/3 (2^/3^). For the functional Bell inequal- 
ity in 14], violation implies A > 2 {2^ /tt^). One can 
easily check that the three Bell inequalities considered 
are inconsistent with the non-bipartite-distillability con- 
dition, Eq. (|19|) . Hence, the violating of these three Bell 
inequalities implies the existence of some bipartite distil- 
lability. 
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